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Abstract 


The  applicability  of  adaptive  array  concepts  to  continuous  aperture 
antennas  was  studied  and  appropriate  aperture  field  distributions  for 
pattern  nulling  were  found  froa  then.  The  adaptive  array  weights  were 
found  to  be  useful  as  discrete  points  in  a  continuous  distribution.  This 
distribution  could  then  be  used  in  an  aperture  integration  scheae  to 
produce  a  nulled  pattern. 

Also  studied  was  the  use  of  geometrical  optics  to  calculate  the 
aperture  field  distribution  of  an  arbitrarily  shaped  reflector.  Under 
soae  restrictions,  geometrical  optics  can  provide  a  useful  approximation. 
Constructing  the  aperture  field  of  a  reflector  defined  by  a  discrete 
grid  of  points  using  a  numerical  ray  tracing  scheae  was  also  investigated. 
Certain  nuaerical  problems  were  identified. 

Finally,  an  attempt  was  aade  to  implement  the  nulled  pattern  by  & 
well  known  beaashaping  method  based  on  geometrical  optics  principles. 

This  technique  was  found  to  be  inadequate.  More  promising  techniques 
for  implementing  the  aperture  distributions  were  suggested  but  not 
pursued  In  this  work. 


I  Introduction 

The  Problem 

Many  systems  which  use  antennas  must  contend  with  the  possibility  of 
jamming  by  hostile  transmitters.  System  designers  often  attempt  to  deal 
with  noise  by  devising  signal  processing  schemes  to  Improve  the  slgnal-to- 
nolse  ratio  of  the  system.  This  may  be  the  only  way  for  most  types  of 
noise,  such  as  thermal  noise  In  the  receiver.  However,  one  often  has  a 
little  more  Information  about  hostile  jammers.  Using  this  information 
presents  the  system  designer  with  yet  another  tool  for  raising  the  signal - 
to-nolse  ratio  of  his  creation. 

Thoughts  like  these  have  led  to  the  development  of  adaptive  antenna 
systems.  Such  an  antenna  modifies  Its  radiation  pattern  according  to 
some  predetermined  criterion  through  feedback  control  (Reference  1).  What 
does  this  accomplish?  Basically,  an  adaptive  antenna  Improves  the  slgnal- 
to-nolse  ratio  of  the  system  It  serves  by  rejecting  Interference.  Some 
adaptive  antennas  use  knowledge  of  the  desired  signal  to  track  that  signal 
and  reject  all  else.  If  Instead,  one  has  knowledge  of  the  spatial  distri¬ 
bution  of  noise  sources  (Ref  1),  he  might  configure  the  adaptive  antenna 
to  place  pattern  nulls  on  those  particular  directions.  Either  way,  one 
seeks  to  maximize  the  signal -to-nolse  ratio  by  rejecting  the  Interference 
with  the  antenna  Itself. 

To  date,  adaptive  antenna  designers  have  confined  their  attention 
primarily  to  array  type  antenna  (Reference  1).  Controlling  the  radiation 
pattern  of  an  array  antenna  Is  easily  accomplished  by  varying  the 
amplitude  and  phase  of  the  Input  signals  which  drive  each  array  element 
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(Reference  l)*  Since  such  control  is  strictly  electronic,  the  system  does 
not  require  mechanical  assistance;  hence,  the  interest  in  array  antennas. 

To  control  the  radiation  pattern  of  a  reflector  antenna  with  a  continuous 
aperture,  one  would  have  to  devise  some  scheme  to  distort  the  reflector 
from  Its  equilibrium  shape  In  a  well  controlled  manner. 

Such  a  scheme  requires  as  a  foundation,  some  knowledge  of  which 
reflector  geometries  are  capable  of  generating  radiation  pattern  nulls 
In  predetermined  locations.  One  must  know  if  such  shapes  exist  and  if 
these  shapes  are  physically  realized.  Clearly,  a  reflector  with  a  loop 
in  it  simply  will  not  do.  This  work  attempts  to  address  the  problem  of 
finding  shapes  that  place  nulls  in  the  radiation  pattern  at  specified 
locations. 

Only  the  simplest  techniques  were  used  In  this  preliminary  work. 
Principles  of  geometrical  optics  were  applied  to  reflector  surfaces  to 
generate  approximate  aperture  field  distributions.  Geometrical  optics 
assumes  that  reflector  dimensions  are  much  greater  than  the  wavelength 
(Ref  12:454  ).  Furthermore,  the  radius  of  curvature  of  the  reflector 
surface  at  any  point  must  also  remain  much  greater  than  the  wavelength 
(Ref  9:  123  ),  This  all  means  that  consideration  here  is  confined  to 
large  reflectors  that  are  not  too  "bumpy".  Once  the  aperture  distribution 
was  determined,  simplified  aperture  integration  generated  the  radiation 
pattern.  These  approximations  mean  that  the  results  are  really  only 
qualitative;  but  they  do  give  some  idea  of  the  location  of  nulls  in  the 
radiation  pattern,  at  least  for  angles  not  far  from  boreslght.  The  analyses 
were  also  limited  to  the  two  dimensional  case,  in  the  Interest  of  simplicity. 

In  the  work  that  follows,  a  brief  description  of  Applebaum  adaptive 
array  antenna  concepts  will  be  presented  in  order  to  familiarize  the 
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reader  with  the  adaptive  array  Ideas  which  are  later  applied  to  apertures. 
In  Section  II,  the  aperture  Integration  technique  used  in  this  work  Is 
described  along  with  proposed  aperture  field  distributions  that  generate 
pattern  nulls.  Section  III  deals  with  geometrical  optics  as  applied  to 
large  reflectors.  The  algorithm  used  to  trace  the  ray  path  from  feed 


to  aperture  Is  presented  also.  Section  IV  attempts  to  implement  a  nulled 
pattern  using  well  known  beam-shaping  techniques  also  based  on  geometrical 
optics.  Section  V  attempts  to  find  reflectors  which  Implement  the  aperture 
distributions  found  In  Section  II.  A  summary  of  the  results  rounds  out 
this  document. 

Possible  Antenna  for  Reflector  Distortion 


At  this  point,  the  reader  might  wonder  about  the  feasibility  of  a 
reflector  that  changes  shape  after  construction  Is  completed.  Jeffrey  H. 
Lang,  et  al  have  proposed  an  antenna  which  might  offer  a  possible  solution 
to  this  problem.  They  have  suggested  a  reflector  constructed  from  a  thin 
membrane  whose  shape  Is  controlled  electronically  (Ref  5,6,7  ). 

Figure  1.1  Illustrates  the  Lang  concept  (Ref  5:  655  )•  An  electron 


beam  controls  the  shape  of  the  reflector  (Ref  5:  655  )•  Referring  to 
Figure  1.1,  the  cocimand  surface  is  pulled  into  a  roughly  parabolic  shape 
by  guy  wires  and  the  reflector  is  pulled  toward  the  command  surface  by 
application  of  a  bias  potential  resulting  In  an  approximately  parabolic 
reflector  (Ref  5:  655).  The  computer  uses  optical  measurements  to 
determine  If  the  reflector  shape  meets  tolerances.  If  not.  It  Instructs 
the  electron  beam  to  scan  the  command  surface  In  such  a  way  that  the  modi 
fled  charge  distribution  and  resulting  electrostatic  force  refine  the 
reflector  shape  (Ref  6:  992).  This  results  In  a  parabola  which  meets 
some  predetermined  tolerance. 
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Fig  1-1.  Lang  Concept 

If  such  an  antenna  were  perfected,  it  might  prove  useful  in  the 
reflector  distortion  problem.  Instead  of  warping  the  membrane  into  a 
parabolic  shape,  perhaps  the  charge  distribution  could  be  altered  to 
form  null  generating  shapes.  Such  a  scheme  would  be  superior  to  purely 
mechanical  means  especially  for  large  reflectors.  J.  H.  Lang,  et  al 
envision  reflectors  up  one  kilometer  in  diameter  (Ref  7  ).  A  mechanical 
scheme  for  such  an  antenna  must  be  clever  indeed.  However,  an  investi¬ 
gation  of  the  electrostatic  reflector's  applicability  to  the  nulling 
problem  Is  beyond  the  scope  of  this  work  and  will  not  be  considered 
further. 

Applebaum  Adaptive  Array 


Applebaum  devised  an  adaptive  array  scheme  which  maximized  the 
single-to-noise  ratio  (Ref  1).  It  accomplishes  this  by  placing  nulls 


in  the  radiation  pattern  on  the  angular  positions  of  any  jamming  transmitter. 
Applebaum  assumed  that  the  angles  of  jammers  was  known  when  he  developed 
his  control  law  for  array  weights. 

Figure  1.2  (Ref  1:  585)  provides  a  starting  point  for  the  description 

of  Applebaum' s  array.  The  array  in  question  has  K  elements.  The  desired 
input  signal  at  element  1  Is  and  the  undesirable  noise  Input  at  element 
1  Is  n.j .  Actually,  these  are  the  complex  envelopes  of  the  desired  signal 
and  noise.  The  output  desired  signal  v  is  the  sum  of  the  x.  weighted  by 
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the  appropriate  w...  The  output  noise  signal  vn  is  the  sum  of  the  n.. 
similarly  weighted  by  the  w.. 

One  knows  from  probability  and  communication  theory  that  p^-  = 

E  [n.*n.]  gives  the  covariance  of  noise  signals  n.  and  n..  In  developing 

I  J  I  J 

his  array  concept,  Applebaum  assumed  that  the  o-fj  are  known.  They  depend 
only  on  receiver  noise  and  the  angular  distribution  of  the  noise  sources 
(Ref  1:  585).  He  then  defined  the  noise  covariance  matrix  as  M  *  [p^ ] 
where  i  and  j  both  vary  from  one  to  K  (Ref  1:  586  ). 

Applebaum  further  assumed  that  the  desired  signals  were  narrowband. 

In  general,  the  interelement  phase  shift  in  an  array  is  frequency 
dependent  (Reference  1).  However,  if  the  desired  signal  has  narrow- 
bandwith,  the  interelement  phase  shift  Is  essentially  constant  and 
depends  only  on  interelement  spacing  and  angle  of  incidence.  This  enabled 
Applebaum  to  write  the  desired  signals  x^  as  x^  *  os^  where  a  is  the  com¬ 
plex  envelope  of  the  desired  signal  Incident  on  the  array  and  the  s^  Is 
a  phase  factor  depending  on  angle  of  incidence,  interelement  spacing,  and 
element  pattern.  Assuming  isotropic  elements,  the  total  phase  shift  at 
the  1^  element  referenced  to  the  first  element  is  BL(i-l)  where  e  *  2-n/x 
and  L  *  d  sin  9  as  shown  in  Figure  1.3.  Therefore,  the  s.  are  given  by 


Fig  1-2.  Linear  Array 


(Ref  1:  586): 


Sj  *  exp  (jedl  sin  e)  ( 

With  these  definitions  in  hand,  Applebaum  determined  that  the 
optimum  weights  which  maximized  signal-to-noise  ratio  are  given  by 
(Ref  1:  586  ): 

MW  -  yS*  ( 


(1-1) 


(1-2) 


where  M  Is  the  noise  covariance  matrix  defined  previously,  W  is  the  vector 
of  array  weights,  S*  is  the  vector  of  S^*  with  give  by  Eq  (  1-1  ). 

The  star  denotes  the  complex  conjugate.  The  y  is  an  arbitrary  constant. 
The  reader  should  refer  to  Appendix  A  for  verification  of  this  result. 


Bore sight 


Applebaum  then  generalized  this  result.  If  no  noise  is  incident, 

then  the  noise  covariance  matrix  reflects  only  receiver  noise  which  is 

equal  in  all  channels  and  uncorrelated  between  channels.  This  defines 

the  quiescent  noise  covariance  matrix  as  M  *  pi.  where  p  is  quiescent 

q  q  k  q 

noise  power  and  Ik  is  a  K  x  K  identity  matrix  (Ref  1:  589).  Since 

Mq"1  *  —  1^,  the  quiescent  weights  are  W  =  vMq_1S  =  S*.  Thus  W  and 

pq  pa  q  Pq  q 

S*  are  scalar  multiples  and  S  =  Wq*  *s  a  measure  of  the  quiescent 
weights  and  hence  quiescent  pattern  (Reference  1).  In  fact,  S  corresponds 
to  uniform  quiescent  weights.  If  one  wants  different  quiescent  weights, 

P 

he  merely  defines  I  s-^  Wq*  with  the  new  weights  and  using  MW  =  MT*  comes 
up  with  the  modified  control  law  (Ref  1:  588  ): 


This  accepts  some  degredation  of  optimal  signal -to-noise  ratio  in  order 
to  achieve  control  over  the  quiescent  pattern  (Ref el:  587  ). 


According  to  Applebaumi  Equation  (1-3  )  gives  weights  that  place  a 
pattern  null  on  a  jammer.  To  illustrate  this,  and  see  what  the  array 
weights  are,  consider  a  K  element  array  whose  mainbeam  is  angle  e  from 
mechanical  boresight.  The  quiescent  noise  covariance  matrix  is 
M  *  p  L  as  stated  before.  The  quiescent  weights  are  given  by: 


& 


Wq  = 


al 

a2  e-jed  sin  es 
a3  e-j2fJdsin0s 


a  e-j(K-l)ed  sin  e$ 
Ic 


(1-4) 


The  a.  are  the  real  number  amplitude  weights  and  if  the  a^  are  all  equal, 
the  quiescent  pattern  is  just  the  familiar  sin(Kx)/sin(x)  (Ref  1:  589  ). 
From  fundamental  antenna  theory,  the  quiescent  pattern  is  (Ref  1:  589  ): 


F  (e)  =  z  a.  e 
q  i=l  1 


j(i-l)8d(sin  e-sin  ©_) 


(1-5) 


where  e  is  the  angle  of  observation  referenced  to  the  mechanical  boresight. 
Vector  notation  greatly  simplifies  the  task.  Therefore,  defining  the 
vector  B  for  the  sake  of  convenience  as  (Ref  1:  589  ): 


1 

gjgd  sin  e 

ej2ed  sin  0 
« 

ej(K-l)Bd  sin  © 


A  W-  . 


.1  a  I  1  i 
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permits  the  quiescent  radiation  pattern  to  be  rewritten  as: 

Fq(Q)  *  BTWq  (1-7) 

With  this  preliminary  out  of  the  way,  one  must  now  seek  an  expression 
for  W.  From  Eq  (1-3)  and  Mq  »  pqIk,  the  W  Is  (Ref  1:  590  ): 

w  •  PqM-1Wq  (1-8) 

One  must  now  find  the  noise  covariance  matrix  and  hence  M“3. 


Define  the  jamming  signal  In  the  first  array  element  as  J(t).  Then 
the  jamming  signal  In  the  1th  element  Is  J^(t)  *  J(t)  exp  [j(1-l)ed  sin  ej] 
The  covariance  between  J^(t)  and  Jk(t)  Is: 


E  CJi*(t)Jk(t)]  -  E[J*(t)e-j(1-1)sd  s1n  sJj(t)e3(K-1)6d  s1n  6j] 

-  E[J*(t)J(t)]  e'3(1‘k)sd  s1n  03 

-  E[|j(t)|2]  e-3(1-k)6d  s1n  03 

=  p .  s1n  ej  (1 

3 


Consequently,  the  covariance  matrix  of  the  janmlng  signal  Is  = 

Pj(exp[-j(1-k)ed  slnQj])  for  1  and  k  from  one  to  K.  is  Hermltlan  with 
all  elements  on  the  same  diagonal  equal.  Its  simple  structure  permits  a 
further  modification  (Ref  1:  590  ).  Defining  the  matrix  (Ref  1:  590  ): 


and  also  defining  U  as  a  KxK  matrix  of  ones  leads  to  (Ref  1:  590  ): 

M.  »  p  H*UH  (1-11) 


Since  the  receiver  and  jammer  noise  are  uncorrelated,  the  total  noise 
covariance  matrix  Is  (Ref  1:  590): 


M  =  M  +  M.  =  pi.  +  p.H*UH  (1-: 

M  J  M  *  J 

Applebaum  gives  M”1  as  (Ref  1:  590  ): 

m’1  -  ~  Cik  -  L  Pj--)  H*UH]  (1-: 

°q  k  pq  +  pj 

Applying  this  expression  to  Eq  (1-8)  leads  to  the  following  equation 
defining  the  weights  (Ref  1:  590  ): 

W  *  W  -  ( - £4“)  H*UHW  (1-: 

q  'Pq  +  Kpj  q 

Some  further  simplification  Is  In  order  (  Ref  1:  590): 


al 

_  oj0d  (S<*n  ©J  ■  sin  0S) 
a2  e 

t  jed  (sin  ej  -  sin  es) 

a3  e 

_j(K-l)  0d  (sin  ej  -  sin  es) 
ak  e 


UHWq  »  Fq(ej) 


(1- 
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H*UHWq  »  Fq(ej) 


(1-17) 


e-j2gd  sin  ej 

e-j(K-l)  gd  sin  ej 

Defining  the  vector  In  Eq  (1-17)  as  BJ  allows  one  to  write  Eq  (1-17)  as 
(Ref  1:  590  ): 

H*UI1Wq  =  Fq(ej)  B*  (1-18) 

Substituting  Eq  (1-18)  into  Eq  (1-14)  results  In  (Ref  1:  590  ): 

w  *  w,  —  Fq<ej)Bj  <!-»> 

To  get  the  nulled  pattern,  one  multiplies  Eq  (1-19)  by  BT  where  B  is 


defined  in  Eq  (1-  6)  (Ref  1:  590  ): 

F(e)  * 

Fq(6)  ' 

H  J 

(1-20) 

BV 

1  +  (sin  0  -  sin  ej)  +  ...  +  gj (K-l) Bd  (sin  e 

-  sin  ej) 

l  ej(K-l)ed  (sin  e  -  sin  ej) 

K=1 

S 

j  [ejed  (sin  e  -  sin  ej)^(K-l) 

(1-21) 

K-l 

Eq  (1-21)  1i 

;  a  geometrical  progression  and  simplifies  to: 

B^B1?  * 

J(K-l)  #  (sin  0  -  sin  ej)  s1n[Kf  (s1n  e  '  s1n  ej)]  (1-2 

D  Dj 

sin[-^-  (sin  e  -  sin 

ej)] 

Since  ( — — )  F_ (©j )  is  just  a  number,  the  second  terra  in  Eq  (1-20) 

Pq  r  *Pj  q 

is  a  sin(Kx)/sin(x)  pattern  with  its  main  beam  at  ej  (Ref  1:  590  ). 
This  comprises  the  nulling  mechanism. 

Eqs  (1-20)  and  (1-22)  provide  an  expression  for  a  pattern  with  a 
null  at  ej.  However,  the  expression  for  the  weights  should  be  broken 
down  for  use  In  Section  II.  Using  Eq  (1-14),  the  1th  weight  Is: 

-  <V,  -  A  l  e'J(,'k)8d  s1n  9j  <V* 

•  (Wq),  -  7^  e-J6d(  s1"  Si  j  [ej6d  sin  ^  ^  (1.2 

Substituting  (Wq)^  *  1  Into  Eq  (1-23)  for  uniform  quiescent  weighting 
gives: 


w  »  !  _  J  e-j6d1  s1n  e  [ e s1n  ej]k 
1  Pq+Kpj  K*1 


Employing  geometrical  progression  and  some  algebra  gives: 


p<  sin  (K  4p  sin  ej)  Bd 

-rfrr-  - £j - J(l+K-21)  2  sin  ej 


Wi  -  1  -  pq+Kp j  sin  (  4r  Sin  ej)  e' 


Eq  (1-25)  will  figure  prominently  in  Section  II  while  Eqs  (1-20)  and 
(1-22)  will  be  Important  in  Section  IV. 
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II  Pattern  Calculations 


Aperture  Integration 


All  the  radiation  patterns  In  this  work  were  computed  by  Integrating 
an  aperture  electric  field  distribution  over  an  aperture  plane.  For 
simplicity*  calculations  Involved  only  two  dimensional  geometries;  there¬ 
fore,  only  one  dimensional,  single  Integrations  were  Involved. 

The  expression  for  radiation  pattern  follows  from  the  following 
definition  for  an  electric  vector  potential: 


E  *  -  v  x  7 


(2-1) 


where  E  Is  electric  field  and  7  is  electric  vector  potential.  The  F 
solves  a  vector  wave  equation  with  a  magnetic  current  source  term 


(  Ref  4:  56  ).  It  is  given  by  (  Ref  4:  56  ): 


(2-2) 


where 

R  -  | r-r1 1 

r  *  Observation  Point  Position  Vector. 
r1  *  Source  Point  Position  Vector. 
M^?1)  ■  Magnetic  Surface  Current. 

z  *  Surface  on  which  1?  (r*)  exists. 


The  far-fleld  approximation  for  this  is  (Ref  12:  379) : 


V 


where  r  Is  the  spherical  coordinate  radial  unit  vector  and  r  is  the  radial 
distance  to  the  observation  point. 


JTOs 

VN* 


In  aperture  integration  problems,  the  magnetic  surface  current  density 
on  surface  z  is  an  equivalent  current  derived  from  the  electric  field  on 
the  surface  (Ref  12:376  ).  The  expression  is  (Ref  12:376  ): 

-  Exn  (2-4) 

A 

where  T  Is  the  electric  field  on  surface  z  and  n  is  the  outward  normal 

* 

unit  vector  to  z.  If  one  selects  a  planar  z,  then  n  remains  constant  on 
the  surface  of  integration.  Substituting  Eq  (2-4)  into  Eq  (2-3)  yields 
fltef  12  :  382  ): 

7  -  |^6r  [  /  l  ej6r,?1  ds1]  x  n  (2-5) 

£ 

The  aperture  surface  z  Is  often  chosen  to  lie  in  the  xy-plane.  The 
aperture  electric  field  Ea  is  also  assumed  to  be  everywhere  tangential 

d 

A 

to  the  aperture  plane  (  Ref  12:381)*  Under  these  conditions,  n  is  simply 
the  z-directed  unit  vector  and  ds*  is  just  dx*dy*  where  the  primes  refer 
to  points  in  the  aperture  plane.  For  convenience,  the  bracketed  term 
In  Eq  (2-5)  will  be  called  P  (  Ref  12:382  ). 

Since  ¥.  has  only  x  and  y  components,  the  P  vector  has  only  x  and  y 

d 

components.  Identifying  these  components  as  Px  and  Py  and  substituting 
into  Eq  (2-5)  results  In  (  Ref  12:382): 

-jer 

7  "  fc:  (py  x  -  Px  y)  (2-6) 

A  A 

where  x  and  y  are  the  x  and  y  directed  unit  vectors  respectively. 

Given  the  magnetic  vector  potential,  defined  by  H  ■  ?  x  J,  the  far 

wv  A  •  • 

zone  electric  field  is  E  *  jw \iK  less  the  r  component,  where  w  is  radian 
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frequency  and  y  Is  the  magnetic  permeability  (  Ref  12:  25).  Since 
T  is  the  dual  of  A,  the  far  zone  H  field  Is  simply  H  *  ju>«F  less  Its 

A 

r  component,  where  6  Is  the  permittivity  (  Ref  12:379).  The  far  zone 
fields  exhibit  plane  wave  behavior;  therefore,  T  ff  x  r  ■ 

-j my/jmT  x  r.  Since  ■  zJ/c»  where  c  Is  the  speed  of  light,  we 
have  *  ^7^  *  8.  This  gives  the  far  zone  electric  field  as 

(Ref  12:  383): 


E  -  jpr  x  F  (2-7) 

Substituting  Eq  (2-6)  Into  Eq  (2-7)  yields  the  expressions  for  far  zone 
electric  field  (Ref  12:  383): 
e-jBr 

E«  *  (Px  cos  0  +  Py  sin  0)  (2-8) 


.  j  01T 

E0  »  jpe4irr  cose  (Py  cos  0  -  Px  sin  0)  (2-9) 

To  develop  the  pattern  equation,  one  must  consider  the  geometry 
Illustrated  in  Figure  2-1.  The  geometry  Is  assumed  to  be  completely 
symmetrical  In  the  x  direction.  The  aperture  plane  lies  on  the  y  axis 
and  extends  from  yamin  to  yamax*  In  this  situation,  0  »  r  I  ^  The  feed 
is  an  electric  line  source  along  the  x  axis  In  the  positive  direction. 

A 

This  results  In  an  x  directed  aperture  field  given  by  Ea  *  x  6(x)  Eax(y). 
Feeding  this  information  into  Eq  (2-9)  gives: 


E0  -  -j£ 


-Jar 
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Eaxty1)  e  s1ne  dy1 


(2-10) 


ram1n 


A  — 1  _  i  -  i 

where  r«r  was  executed  using  r  ■  yy  . 


FEED i  X -DIRECTED  ELECTRIC 
LINE  SOURCE 


1“  y«*in 


FLg  2-1.  Apert  ur« 


When  pattern  calculations  are  made,  the  radial  distance  remains 
constant  and  terms  Involving  It  are  eliminated  along  with  other  constant 
terms  by  normalization.  The  resulting  for  field  pattern  Is  given  by: 


G(e) 


'amax 


cose 


+  ttllAJA  4 

/„  Eax(yl)  s1ne  dy1 


ram1n 


(2-11) 


'amax  , 

maxlcose/  Eaxfy1)  eJ6y  sine  dy1 
yam1n 


When  using  aperture  Integration,  one  normally  confines  his  attention  to 
angles  near  boreslght  where  cost  will  remain  close  to  unity.  Conse¬ 
quently,  the  final  expression  for  pattern  used  In  this  work  is: 


(2-12) 


Figure  2-2  shows  Eq  (2-12)  applied  to  a  uniform  distribution  over  a 
100  wavelength  aperture.  Figure  2-3  shows  the  results  from  a  10  wave¬ 
length  aperture.  Analytic  integration  results  in: 


6(e)  - 


sin 


(  sin  e 


■y-  sin  e 


(2-13) 


For  the  10  wavelength  aperture,  the  first  null  occurs  at  sin  0  ■  0.1 
corresponding  to  e  *  5.7  degrees. 

In  the  preceding  discussion,  the  aperture  fields  were  assumed 
tangential  to  the  aperture  plane,  that  Is  the  field  vectors  lie  completely 
In  the  plane.  For  non-uniform  phase  distributions,  however,  this  is  not 
the  case  (Ref  10:  162).  The  direction  of  propogation  of  the  reflected 
wavefront  Is  normal  to  the  elkonal  surfaces.  Since  in  the  geometrical 
optics  approximation  the  fields  exhibit  plane  wave  behavior,  the  H  field 
will  not  lie  In  the  aperture  plane  for  the  x  directed  E  field  considered 
previously  If  the  phase  distribution  varies  over  the  aperture. 

If  one  considers  the  effect  of  non-uniform  phase  distributions,  he 
must  employ  an  alternative  expression  such  as  In  Reference  ten.  The 
expression  for  far-zone  electric  field  is  (Ref  10:  161): 

T  *  r  x  J {n  x  Ta  -  [r  •  (s  x  la)  n  -  (s  x  TTa)(n  •  r)]}e^6r'r  ds1 

z 


(2-14) 


gp  oC- 


where 


z  =  plane  aperture  surface 

A 

n  =  outward  normal  to  z 

A 

r  =  unit  vector  pointing  from  origin  to  field  point 
r*  *  vector  from  origin  to  ds*  in  aperture 

A 

s  a  unit  vector  along  ray  trough  aperture 
Eq  (2-14)  reduces  to  Eq  (2-12)  if  the  two  dimensional  assumptions  are 

A 

applied  and  the  phase  deviation  from  uniform  are  assumed  small  giving  s 

A 

approximately  equal  to  n  (Ref  10:  162).  Ignoring  the  phase  error  is  not 
really  correct  but  is  frequently  done  to  simplify  calculations;  however, 
the  results  have  primarily  qualitative  value  (Ref  10:  173). 

Applebaum  Weights  and  Patterns 

One  reason  for  using  the  simplified  aperture  integration  formula  is 
its  similarity  to  the  summation  expression  found  in  array  calculations. 
This  suggests  that  one  might  use  the  Applebaum  array  weights  as  discrete 
points  in  an  aperture  distribution.  Since  the  aperture  integral  is 
analogous  to  the  array  summation,  such  a  distribution  should  produce  a 
similarly  nulled  pattern.  This  is  in  fact  the  case. 

To  illustrate  this  Eq  0.-25)  from  Section  I  was  applied  with  K  *  201 
and  d  *  0.05  on  a  ten  wavelength  aperture.  Figure  2-4  shows  the  aperture 
magnitude  distribution  and  Figure  2-5  shows  the  phase  distribution  for 
creating  a  null  in  Figure  2-3  at  two  degrees.  Similarly  Figures  2-6  and 
2-7  show  aperture  magnitude  and  phase  respectively  for  a  null  at  seven 
degrees,  while  figures  2-8  and  2-9  correspond  to  a  null  at  ten  degrees. 
Figures  2-10  through  2-12  show  the  patterns  resulting  from  the  appropriate 
distributions  for  two,  seven,  and  ten  degrees  respectively. 
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Further  experimentation  indicates  that  the  phase  distribution  is  the 
dominant  factor  in  determining  null  location.  Figure  2-13  shows  that 
aperture  integrating  the  phase  distribution  of  Figure  2-5  with  a  uniform 
amplitude  distribution  creates  a  null  at  2.7  degrees.  Similarly,  Figures 
2-7  and  2-9  integrated  with  uniform  amplitude  distributions  yield  nulls 
at  6.7  degrees  and  10.7  degrees  as  shown  in  Figures  2-14  and  2-15. 
Implementing  the  amplitude  distributions  with  flat  phase  does  not  con¬ 
sistently  come  as  close.  Such  calculations  result  in  nulls  at  4.5,  6.3, 
and  10.9  degrees  for  the  three  respective  cases.  Especially  in  the  two 
degree  null  case,  the  flat  phase  patterns  do  not  deviate  much  from  the 
quiescent  pattern  with  nulls  at  5.7  and  11.5  degrees  as  shown  in  Figures 
2-16,  2-17,  and  2-18.  This  implies  that  the  phase  distribution  should 
receive  priority  in  any  attempt  to  implement  the  aperture  distributions 
shown  in  Figures  2-4  to  2-9.  Using  only  the  phase  distribution,  however, 
results  in  an  as  yet  undetermined  error  which  deserves  further  study. 


Ill  Secondary  Pattern  Calculation 


Geometrical  Optics 

In  order  to  approximate  the  field  distribution  in  the  aperture, 
geometrical  optics  principles  were  applied  to  the  reflector  shapes  in 
question.  The  use  of  geometrical  optics  requires  that  the  dimensions 
of  the  scattering  surface  be  large  compared  to  a  wavelength  (Ref  12:454  ). 
It  also  places  a  constraint  on  the  smoothness  of  the  surface,  since  the 
radius  of  curvature  of  the  surface  must  always  be  large  compared  to  a 
wavelength  (Ref  9:  123  ).  Even  if  these  requirements  are  met,  however, 
the  approximation  is  still  rather  crude.  All  diffraction  effects  have 
been  ignored.  Thus,  the  reader  should  interpret  the  results  qualitatively. 
The  null  positions  and  the  effect  of  reflector  shape  on  them  are  the  points 
to  watch. 

Specialized  to  the  two-dimensional  case,  the  geometrical  optics 
expression  for  amplitude  is  (Ref  12:  450): 

|E| '  |E°'  V"v^  (3'1) 

where  |Eq|  is  the  magnitude  of  the  field  at  the  reference  location  and 
r  is  the  radius  of  curvature  of  the  wavefront  at  the  reference  point. 

\*TI 

The  p  denotes  the  distance  from  the  reference  point  to  the  field  point. 

To  Include  phase,  one  simply  applies  a  phase  factor  determined  by  p. 

The  resultant  expression  is  (Ref  12:451  ): 


where  e  is  simply  2ir/x  and  0Q  is  the  reference  phase  at  the  reference 
point. 

In  order  to  compute  the  aperture  field  contribution  of  a  single  ray, 
one  must  first  determine  the  field  radiated  from  the  feed  at  the  point  of 
impact  for  that  ray  on  the  reflector.  In  this  work,  the  feed  is  assumed 
to  be  a  line  source  at  the  origin.  Placing  the  reference  point  one  wave¬ 
length  from  the  feed  and  choosing  |EQ|  =  1,  0Q  =  -e  in  Eq  3-2  yields: 


I  1 

1  +  p 


,-jed  +  p) 


(3-4) 


for  the  field  phasor  at  the  point  of  reflection.  But  1  +  p  is  simply 
the  distance  from  feed  to  reflection  point,  so  letting  1  +  p  *  rj  gives: 


V-7 


e-j6rl 


(3-5) 


To  calculate  the  field  phasor  in  the  aperture,  one  must  compute  2  ,  the 
distance  from  the  reflection  point  to  the  point  where  the  ray  in  question 
intersects  the  aperture.  Then  using  |EQ|  =  and  0O  "  "Brl  3"^ 
gives: 


!  cw 
r„.  +  r. 


,-jer2  g-jBrj 


(3-6) 


All  that  remains  in  finding  rcw>  the  radius  of  curvature  of  the  reflected 
wavefront.  This  is  given  by  (Ref  12:452  ): 


_JL  =  —  + 

rcw  rl 


r  cos 
c  o 


(3-7) 


where  rc  is  the  physical  radius  of  curvature  of  the  reflector  at  the 
reflection  point  and  0Q  is  the  angle  of  incidence  of  the  incoming  ray 
with  respect  to  the  surface  normal.  If  the  surface  is  convex  when 


viewed  from  the  line  source  feed,  r  is  greater  than  zero.  If  the 
surface  is  concave,  the  r  is  less  than  zero.  Thus,  given  a  ray,  Eq  3-6 
determines  the  aperture  field  due  to  it. 

Ray  Tracing 

Before  one  can  compute  the  aperture  field  approximation,  he  must 
first  trace  given  rays  from  the  feed  through  their  points  of  reflection 
to  their  aperture  intersection  points.  Figure  3-1  illustrates  the 
geometry.  The  e  is  the  angle  the  incident  ray  makes  with  the  horizontal 
axis.  The  reflector  surface  is  described  by  a  discrete  grid  of  equi- 
spaced  points. 

The  first  step  is  to  find  the  (x^,  y..)  impact  point  from  ©^  and  the 
grid.  The  line  y  *  kx  represents  the  incident  ray  in  this  two  dimensional 
geometry,  where  k  *  tan©^.  To  determine  the  impact  point,  one  must  define 
a  miss  function  (Ref  11:1449)  gn  =  k*n-yn*  Testing  each  grid  point  will 
yield  same  yn-1)  and  (*nlyn)  between  which  the  miss  function 

changes  sign,-  because  9n_i9n  <  0.  The  impact  point  lies  on  the  reflector 
surface  between  these  two  points.  This  impact  point  is  determined  by 
(Ref  11:1449): 


X1  *  yn(1-F)  +  Xn-1F 

yi  '  Vl-F)  +  VlF 


(3-8) 


where  F  -  gn  /(g^-g^!). 

To  determine  the  point  of  impact  in  the  aperture  plane,  one  must 
consider  Figure  3-2.  The  0  and  the  frare  the  angles  made  by  the  surface 
normal  and  the  reflected  ray  respectively  with  the  horizontal  axis. 

Snell's  Law  at  the  impact  point,  valid  for  the  smooth  surfaces  in  question. 


Fig  3-1.  Incident  Rays 

dictates  that  the  angle  of  incidence  and  angle  of  reflection  with  respect 
to  the  surface  normal  be  the  same  (Ref  11:1449).  This  angle  is  denoted 
by  <5  in  Figure  3-2.  Simple  geometry  then  yields: 

x  *  20  -  9  (3-9) 

Determining  0  requires  approximation  of  the  first  derivative  of  the 
reflector  surface  at  the  impact  point.  Calculating  it  at  (xnl ,yn_1) 
and  (xn»yn)  ^  a  we^  established  finite  difference  formula,  then  using 

linear  interpolation  to  (x.j,y^)  yields  acceptable  results.  Once  this 

derivative  Is  determined,  0  can  be  found.  The  tangent  unit  vector  to 

A  A  A 

the  surface  at  (x.  is  given  by  t  *  cosax  +  sin  ay,  where  a  is  the 
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Fig  3-2*  Angle  Definitions 


angle  made  by  t  and  the  horizontal.  The  slope  of  the  line  determined 
by  t  Is  tana  and  since  this  line  is  the  tangent  line  to  the  curve  at 
(x1,yi): 

*"»  *  $  I  X  .  X,  <3*10> 

where  dy/dxjx=x..  was  already  discussed.  The  normal  vector  at  (x-.y^)  is 

*  A  A 

n  *  dt/da  *  sinax  +  cosay.  The  slope  of  the  line  determined  by  the 
normal  is  tan0;  hence,  tan0  *  cosa/(-sina)  =  -cota  which  implies  that 
tana  *  -cot0,  Substituting  back  into  Eq  3-20  yields  0  in  terms  of  a  known 
quantity  (Ref  11:1449  )i 

-c°t0  *  x  *  x^  (3-11) 
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One  must  consider  some  computational  points  before  actually  using 
Eq  3-11  to  find  0.  The  range  of  0  is  limited  to  0  £  0  <_  it.  If 
dy/dx|x=x.j  >  0,  then  0  >  * / 2*  If  dy/dx|x=x.j  <  0,  then  0  <  The 

case  dy/dx|x*Xj  *  0  implies  that  0  =  7r/2*  Any  scheme  to  compute  0 
should  reflect  this.  Since  arctangent  implementations  on  computers 
generally  range  only  from  -Vg  to  one  must  exercise  care  in  the 
computation.  A  suitable  scheme  is  Eq  (3-12). 


t*"'1n$ix.x1>'1n  £ix.Xl<° 


»■(-?- 


*  -  ^-1[(flx.xrl] 


dy 
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(3-12) 


Using  0  in  Eq  3-9  gives  From  Figure  3-3,  the  point  of  intersection 
of  the  reflected  ray  with  the  aperture  (x  ,0)  is  given  by: 


xi  -  yi/tan  ft 


x  2 


X  2 


(3-13) 


V  xi  +  y^/tan  (w-ffr)  x  >  | 

To  compute  the  field  contribution,  one  simply  applies  Eq  3-6  along 
with  Eq  3-7.  The  physical  radius  of  curvature  at  (x^,y^)  of  the  reflector 
surface  is  simply: 


O  *  Ix.x/J3/2 

dx  x*Xj 


(3-14) 
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where  the  sign  of  the  second  derivative  is  preserved  in  order  to  make  r£ 
appropriately  positive  or  negative  based  on  convexity  or  concavity. 

The  process  of  approximating  the  second  derivative  exhibits  an 
inherent  Instability  which  defied  efforts  at  eliminating  it.  The  effect 
of  the  Instability  will  show  up  graphically  in  the  spherical  reflector 
example.  Appendix  B  describes  the  problem. 

Aperture  Field 

After  defining  an  aperture,  usually  centered  on  the  reflector's 
boresight  and  wide  enough  to  include  all  reflected  ray  Intersection 
points,  one  selects  a  set  of  discrete  points  in  the  aperture  at  which 
he  wishes  to  compute  the  field.  Taking  each  point  in  turn,  one  determines 
which  rays  Intersect  the  aperture  at  that  point  and  adds  the  field 


contributions  of  these  rays.  For  most  relatively  smooth  reflectors,  only 
one  ray  will  strike  the  point  in  question. 


To  find  the  right  ray,  one  steps  through  possible  values  of  0  until 
he  finds  e,  and  e0  such  that  x„0  <  x.  <  x.  where  x.  is  the  desired 
aperture  point.  Then  one  computes  x^  from  ©3  =  (©j  +  ©2)/2.  Then  one 
uses  a  bisection  algorithm  to  close  in  on  the  correct  value  of  0.  If 
xfl  -  xfl3  <  0,  then  the  correct  0  is  between  ©3  and  ©2.  If  xfl  -  xfl3  >  0, 
then  the  0  is  between  ©j  and  ©3-  The  bisection  process  is  then  repeated. 
When  xa  -  xa3  is  sufficiently  close  to  zero,  the  process  terminates. 

When  completed,  the  process  yields  a  discrete  grid  of  aperture  field 
points. 

Parabolic  Example 

The  scheme  described  before  will  now  be  tested  on  a  parabola. 

Figure  3-4  shows  the  parabola.  Its  diameter  is  100X  where  x  is  the 
wavelength.  The  focal  length  is  50x  with  focus  at  the  origin.  The  feed 

A 

is  an  x  directed  electric  line  source  at  the  origin.  The  aperture 

A 

electric  field  will  also  be  x  directed,  by  application  of  Snells  Law 
(  Ref  12:426  ). 

Analyzing  the  parabola,  one  notes  that  the  phase  will  be  constant 
on  an  aperture  centered  on  boresight  and  100X  wide.  Thus,  the  wave  front 
Incident  on  this  aperture  is  a  plane  wave.  Since  a  plane  wave  has  constant 
amplitude,  attenuation  occurs  only  between  the  feed  and  the  reflector. 

Using  this  in  Eq  3-6  yields: 


A  little  analytic  geometry  shows  that  =  50[1  +  (jot))2]-  Normalizing 
this  to  a  unity  maximum  gives: 


1 


Figure  3-5  shows  the  amplitude  distribution  generated  by  the  method  of 
this  section.  The  method  duplicates  the  analytic  geometrical  optics 
analysis.  It  also  duplicates  the  constant  phase,  but  the  flat  line 
has  not  been  plotted.  Figure  3-6  shows  the  Section  II  pattern  result. 

Figures  3-7  through  3-9  correspond  to  Figures  3-4  through  3-6, 
but  with  diameter  10\  and  focal  length  5x.  For  such  a  reflector,  the 
validity  of  geometrical  optics  is  questionable,  since  diffraction  effects 
are  present.  However,  for  angles  near  boresight,  the  pattern  results 
are  still  useful. 

Spherical  Example 

Figure  3-10  shows  the  test  surface.  The  radius  of  curvature  is 
200x.  The  diameter  of  the  reflector  is  100a.  The  feed  is  situated  on  the 
boresight  50X  from  the  center  of  the  reflector  and  coincides  with  the 
origin.  The  aperture  is  centered  on  a  feed  and  is  lOOx  wide  to  include  all 
rays. 

Figures  3-11  and  3-12  show  the  amplitude  and  phase  calculated  by 
the  methods  of  this  section  with  analytic  determination  of  the  second 
derivative.  Figures  3-13  and  3-14  shows  the  results  using  the  final 
second  derivative  approximation  in  Appendix  B.  While  the  phase  result  is 
good,  the  amplitude  exhibits  a  messy  error  distribution. 


The  net  effect  of  the  error  is  not  overly  drastic.  The  messy  ampli¬ 
tude  distribution  seems  to  cluster  about  the  correct  one.  Also,  phase  is 
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dominant  in  determining  null  location.  A  comparison  of  patterns  is  in 
order.  Figure  3-15  shows  the  pattern  resulting  from  the  analytic 
derivative  distribution.  Figure  3-16  shows  the  pattern  of  the  messy 
distribution.  The  amplitude  error  has  effected  side-lobe  levels  and 
null  depths  but  not  null  locations.  The  side-lobe  level  and  null  depth 
errors  are  noticeable,  but  not  intolerable. 

Figures  3-17  through  3-23  we  correspond  to  Figures  3-10  through  3-16 
for  a  spherical  reflector  with  all  the  previously  described  dimensions 
reduced  by  a  factor  of  ten.  Comparing  Figures  3-22  and  3-23  shows  that 
null  location  has  not  been  altered.  Again,  the  effect  of  the  amplitude 
error  is  limited  to  side-lobe  level  and  null  depth.  Consequently,  this 
work  will  continue  to  employ  the  technique  of  this  section,  even  though 
the  results  are  only  qualitative,  except  for  null  location.  Clearly, 
future  work  must  refine  the  second  derivative  approximation. 


IV  Beamshapinc 


Reflector  Equation 

The  simplest  method  for  implementing  a  given  antenna  pattern  applies 
geometrical  optics  to  a  known  power  pattern  to  get  a  reflector  shape  that 
generates  the  pattern  (Ref  10:  497).  The  technique  is  crude,  but  might 
provide  a  first  approximation  (Ref  10:  497).  If  the  reflector  resulting 
from  the  method  is  physically  unreal izeable,  one  cannot  however,  simply 
call  the  pattern  unreal izeable.  The  method  is  not  that  conclusive.  The 
reader  should  also  note  that,  because  the  method  is  based  on  geometrical 
optics,  only  reflectors  large  compared  to  a  wavelength  should  be  considered. 

Assuming  that  the  desired  power  pattern  is  P(y)  between  y ^  and  y2» 
one  can  generate  a  reflector  as  follows.  Referring  to  figure  4-1  (Ref  10: 
498),  positive  angles  are  clockwise  from  the  horizontal  and  the  distance 
from  the  feed  point  to  the  reflector  surface  for  a  given  value  of  x  is 
defined  to  be  r(x).  The  x  is  the  angle  from  horizontal  to  the  incident 
ray,  while  y  is  the  angle  from  horizontal  to  the  reflected  ray.  The  y-^ 
and  $2  are  the  limits  of  the  range  over  which  one  wished  to  specify  P(y); 
the  pattern  outside  the  range  will  be  ignored,  but  the  reader  should  be 
aware  that  diffraction  effects  will  yield  non-zero  field.  The  pattern 
level  will  fall  quickly  outside  the  specified  range  as  long  as  the  size 
of  the  reflector  meets  geometrical  optics  criteria  (Ref  10:  497).  The  Xj 
and  X£  define  the  angular  extent  of  the  reflector.  The  geometry  of  figure 
4-1,  coupled  with  the  law  of  reflection,  leads  to  the  equation  defining 
the  reflector  contour  (Ref  10:  498): 


1  dr  x  -  y 

T  dJ  =  tan— T“ 


(4-1) 
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Before  one  can  integrate  Eq  (4-1),  he  must  first  find  the  relation¬ 
ship  between  y  and  x.  This  arises  from  power  considerations.  Geometrical 
optics  demands  that  the  power  in  a  given  tube  of  rays  must  remain  constant 


(Ref  2:  404).  Consequently,  if  no  loss  occurs  at  reflection,  the  power 
in  the  incident  ray  tube  must  equal  that  in  the  reflected  ray  tube  (Ref  10: 
499).  Defining  the  incident  power  pattern  from  the  feed  as  I(x)  and  con¬ 
sidering  incident  and  reflected  ray  tubes  with  angular  width  dx  and  dy 
respectively,  energy  conservation  gives  (Ref  10:  499): 

I(x)  dx  =  K  P(y)  dy  (4-2) 

where  the  K  is  a  constant  of  proportional ity.  The  fact  that  the  total 
power  from  the  feed  must  equal  the  total  reflected  power  determines  the  K 
from  (Ref  10:  499): 


x2 

J  I(u)  du  -  K 
X1 

To  determine  the  relationship  between  y  and  x  from  Eq  (4-2),  one  must 
consider  that  the  power  in  any  feed  range  x^  to  x  will  equal  the  power  in 
some  corresponding  reflected  range  yj  to  y.  Applying  this  information  to 
Eq  (4-2)  yields  (Ref  10:  499): 

y  1  r  X 

P(u)  du  *  —  /  I(v)  dv  (4-4) 

yi  xj 


r 


P(v)  dv 


(4-3) 


Eq  (4-4)  defines  y  implicitly  as  a  function  of  x.  In  order  to  simplify  the 
calculations,  I(x)  =  1  was  assumed,  implying  an  isotropic  feed  over  the 
range  x^  to  Xg.  Using  this  in  Eqs  (4-3)  and  (4-4)  gives  a  simpler  expression 
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for  the  right  hand  side  of  Eq  (4-4): 

/  p(u)  du  =  X  "  Xl  f  2  P(v)  dv  (4-5) 

y,  x2-  xj  J 

n 

One  can  think  of  Eq  (4-5)  as  defining  g(y)  =  f(x).  The  y  can  be  computed 
from  this  as  y  =  g~*(f(x))  which  is  y  as  a  function  of  x.  In  general,  y 
must  be  found  through  interpolation  of  numerical  data,  since  the  analytic 
integration  of  P(y)  may  not  always  be  feasible. 

To  illustrate  the  use  of  this  relationship,  a  simple  example  might 
prove  useful.  If  one  defines  P(y)  =  csc2(y)  cos(y)  then  is  follows  that: 


r 


P(u)  du  = 


sin(y:) 


sin(y2) 


*  pr  (4-6) 


Similarly: 


/ 


P(u)  du  = 


sin(y,) 


sin(y) 


(4-7) 


'y1  '"/li 

Sticking  Eqs(4-6)  and  (4-7)  into  Eq  (4-5)  yields  an  expression  for  y: 


sin 


-1 


(— — P 

\sin(y1)  r  \x2-  xj J  _ 


(4-8) 


This  can  be  used  in  solving  for  r(x). 

Eq  (4-4)  enables  one  to  integrate  Eq  (4-1)  directly.  Doing  so  results 
in  an  expression  for  the  radial  distance  from  the  feed  to  the  reflector 
surface  for  a  given  angle  x  (Ref  10:  499): 
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Fig  4-1.  Geometry  for  Beamshaping 


ln  ~  ■  /Xten[^l  du 

»  • 


o  0 


(4-9) 


This  leads  directly  to: 


r(x)  =  rQ  exp  i 


tan 


u  -  y(u) 


du 


(4-10) 


The  rQ  is  a  scale  factor  which  determines  the  absolute  dimensions  of  the 
reflector.  In  choosing  rQ,  one  must  always  remember  that  the  contour 
was  derived  from  geometrical  optics  principles  and  one's  confidence  in  it 
should  deteriorate  as  the  size  of  the  reflector  decreases.  Also,  the 
technique  Ignores  any  phase  considerations  and  is  crude.  It  is  really  best 
applied  to  attempts  to  Implement  main-beam  shapes,  not  patterns. 


Implementation 

As  stated  earlier,  an  isotropic  feed  simplifies  calculations,  so  this 
was  assumed.  Therefore,  Eq  (4-5)  applies.  Furthermore,  and  x2  were 
chosen  as  -60  degrees  and  +60  degrees  respectively.  This  gives  x2  -  x^  120° 
and  (x  -  x1)/(x2  -  Xj)  =  (3x  +  TT )/(2rr ) .  Given  this,  the  computerized 
application  of  beamshaping  is  simple: 


1.  Select  y2,  y2  ,  and  P(y). 

2.  Compute  Pr  =  y  2P(u)  du  . 

3.  Define  f(x)  =  Pr  (3x  +7T)/(27T). 

y 

4.  Tabulate  g(y)  =  f  P(u)  du  for  y  to  y^. 

y i 

5.  Compute  y  =  g“^(f(x))  by  interpolation  whenever  needed. 

6.  Tabulate  r(x)/rQ  from  Eq  (4-10). 


7.  Compute  rQ  from  chosen  reflector  diameter  and  choice  for  or 

8.  Scale  r(x)/rQ  to  get  r(x). 


To  illustrate  beamshaping,  a  sin(x)/x  example  was  computed.  For  this 
case,  the  expression  for  P(y)  was: 


P(y)  = 


nl  -i  2 

si n(/J2  sin  y) 


sin  y 


(4-11) 


where  L  was  the  diameter  of  the  desired  reflector,  chosen  to  be  100  wave¬ 
lengths.  Figure  4-2  shows  the  result  which  is  discernably  parabolic. 
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Pattern  Null  Attempt 

In  order  to  find  a  nulled  P(y)  to  use  in  the  method  of  this  section, 
the  reader  should  refer  to  Eqs  (1-20)  and  (1-22).  These  equations  define 
the  nulled  pattern  resulting  from  an  Applebaum  array  of  K  elements.  With 
uniform  quiescent  weighting,  the  quiescent  pattern  is  a  sin(Kx)/sin(x) . 

For  a  continuous  aperture  with  a  uniform  distribution  on  it,  the  pattern 
is  a  sin(x)/x.  One  might  form  a  loose  connection  between  the  two  cases 
and  use  the  Sa(x)  =  sin(x)/x  in  place  of  the  sin(Kx)/sin(x)  in  equation 
(1-20)  when  trying  to  find  a  nulled  pattern  for  the  case  of  apertures. 
Interpreting  Kd  as  L  and  using  Sa(x)  instead  of  sin(Kx)/sin(x)  in  Eqs  (1-20) 
and  (1-22)  yields: 

P(y)  =  | Sa(0)  -  exp(jz)  Sa(0j)  Sa(z)  |2  (4-12) 

where 

0  sin  y 

S1"  *s 

z  =Pfysin  y  - sin  yj) 

yj=  angular  location  of  desired  null 
The  reflector  resulting  from  this  is  shown  in  figure  4-3.  The  reader  will 
note  that  the  shape  is  the  same  as  that  for  the  sin(x)/x  pattern  shown  in 
figure  4-2. 

The  beamshaping  technique  cannot  adequately  cover  the  nulling  case, 
because  it  ignores  all  phase  considerations.  Its  primary  application  is  in 
the  synthesis  of  reflector  shapes  which  will  produce  a  given  main-beam 
configuration,  such  as  the  csc-squared  beam  discussed  earlier.  Pursuing 
it  further  is  not  recommended. 
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V  Suggested  Extensions  and  Summary 

Refined  Analysis 

Before  one  can  refine  the  analysis  technique  given  in  this  work,  he 

must  first  solve,  or  at  least  circumvent,  the  second  derivative  instability 

described  in  Section  III  and  Appendix  B.  A  possible  fix  might  be  to  alter 

the  computational  perspective  used  in  the  program  in  Appendix  D.  Instead 

of  calculating  r  as  an  intermediate  step,  one  might  do  better  by  finding 
c 

r£w  in  one  step: 


where  r-  and  rcw  are  the  physical  radius  of  curvature  of  the  reflector  at 
x^  and  the  radius  of  curvature  of  the  reflected  wavefront  respectively. 

The  rj  is  the  distance  from  the  feed  to  the  reflection  point  (x^ ,y^ )  and 
d  is  the  angle  of  reflection  with  respect  to  the  surface  normal.  Such  a 
computational  trick  might  smooth  the  results.  The  idea  is  due  to  Capt. 
Thomas  Johnson,  who  was  teaching  at  the  Air  Force  Institute  of  Technology 
at  the  time  of  this  writing. 

Given  that  one  finds  a  workable  fix,  he  might  seek  an  analysis  tech¬ 
nique  which  more  faithfully  reflects  the  actual  fields  and  patterns.  Much 
of  the  work  discussed  previously  involved  reflectors  ten  wavelengths  across. 
One  can  justifiable  question  the  accuracy  of  geometrical  optics  alone  in 
such  a  case,  since  diffraction  effects  should  be  quite  noticeable.  Fig  3-18 
shows  a  discontinuous  aperture  field  which  is  common  in  a  geometrical 
optics  analysis.  At  ten  wavelengths,  however,  one  expects  the  field  to 
vary  much  more  smoothly.  Generating  the  aperture  field  distribution  with 


First  Diffraction  Point 


Fig  5-1 •  Simple  GTD  Geometry 

a  simple  geometrical  theory  of  diffraction  (GTD)  approach  will  do  much  for 
this  problem.  Such  a  technique  could  also  be  used  to  compute  the  far-field 
pattern;  however,  a  GTD  generated  aperture  field  coupled  with  aperture 
integration  would  be  more  instructive  in  the  nulling  problem. 

Figure  5-1  shows  the  simplest  GTD  model  for  the  problem  under  consider 
ation  here.  The  total  field  at  a  point  on  the  aperture  is  the  sum  of  the 
contributions  from  all  reflected  rays  and  the  two  rays  singly-diffracted 
from  the  edges  (Ref  12:  484).  The  reflected  rays  are  found  by  the  methods 
of  Section  III.  The  calculation  of  the  diffracted  fields  is  discussed 
elsewhere  (Ref  12:  458-477);  a  straightforward  method  specialized  to  the 
problem  at  hand  can  be  found  in  Appendix  C  for  the  reader's  convenience. 
Higher  order  effects  decrease  in  significance  with  increasing  reflector 
size,  inclusion  of  these  is  a  discretionary  matter. 


Suggestion  for  Implementing  Distributions 

The  most  promising  technique  for  implementing  the  distributions  of 
Section  II  is  due  to  Galindo  (Ref  2).  The  method  employs  the  principles 
of  geometrical  optics  to  implement  arbitrarily  specified  aperture  amplitude 
and  phase  distributions  with  a  dual  reflector  geometry  as  shown  in  figure 
5-2  (Ref  2:  404).  The  problem  eventually  reduces  to  the  numerical  solution 
of  a  coupled  pair  of  first  order,  non-linear,  ordinary  differential  equations 
(Ref  2:  404).  The  development  is  intricate  and  the  interested  reader 
should  look  to  Reference  2  for  an  explanation  of  the  method.  This  method 
is  designed  for  dual  reflectors  like  that  in  figure  5-2;  however,  one 
might  be  able  to  apply  the  principles  behind  the  analysis  to  a  single 
reflector  and  develop  a  similar  method,  one  that  considers  phase.  Reference 


8  presents  a  method  based  on  Galindo's  but  computationally  simpler;  however. 


its  application  to  single  reflectors  does  not  look  too  promising. 


Summary 

Aperture  integration  provides  a  viable  technique  for  considering 
pattern  nulling  with  reflector  antennas.  It  is  directly  analogous  to 
the  summation  of  element  field  contributions  found  in  antenna  array 
calculations.  Consequently,  the  array  weights  given  by  Applebaum's 
adaptive  array  formulation  can  be  used  as  discrete  points  in  a  continuous 
aperture  distribution.  Section  II  illustrated  how  these  distributions 
resulted  in  nulled  patterns. 

In  analyzing  large  reflectors  with  non-traditional  contours,  geometrical 
optics  is  quite  useful.  The  usefulness  of  the  method  presented  in  Section 
III,  however,  is  limited  by  the  numerical  instability  of  the  second  deriv¬ 
ative  approximation.  The  reader  will  note  that  the  method  performed  better 
on  the  ten  wavelength  case  than  on  the  one  hundred  wavelength  cases  from 
a  numerical  standpoint.  This  improvement  results  from  using  the  same 
number  of  discrete  points  in  the  ten  wavelength  cases  as  were  used  in  the 

hundred  wavelength  cases,  giving  a  smaller  step  size  for  the  reflector 

grid.  Unfortunately,  as  in  most  numerical  analysis  problems,  further 
reducing  the  step  size  does  not  always  give  better  results.  Increasing 
the  number  of  grid  points  in  the  ten  wavelength  examples  in  Section  III 

in  order  to  shrink  the  step  size  actually  makes  the  problem  worse.  The 

fix  given  at  the  beginning  of  this  section  seems  the  best  candidate  for 
a  solution.  Given  that  the  numerical  problem  can  be  removed,  the  ray- 
tracing/geomstrical  optics  analysis  could  be  useful  in  further  work. 

Coupling  it  with  a  formulation  of  geometrical  theory  of  diffraction  can 
make  it  even  more  useful. 

The  beamshaping  method  presented  in  Section  IV  is  already  well  known. 

It  is  easy  to  implement,  but  too  simplistic  for  the  problem  under  consideration. 


The  results  are  unacceptable  when  applied  to  pattern  nulling,  possibly 
because  all  phase  considerations  are  ignored.  The  method  was  originally 
designed  to  synthesize  main-beam  shapes  and  is  out  of  its  depth  in  the 
nulling  problem.  Galindo's  method  should  provide  a  good  baseline  for 
further  studies.  Other  methods  for  dual  off-set  reflectors,  also  based 
on  geometrical  optics,  may  also  prove  useful. 


Appendix  A:  Applebaum  Control  Equation  Verification 


Verifying  the  Applebaum  control  equation  demands  that  one  first  find 
an  expression  for  the  signal-to-noise  ratio  and  then  show  that  the  control 
law  weights  in  fact  maximize  the  S/N.  As  a  first  step,  one  computes  the 
outputs  vs  and  vn  (Ref  1:  586  ): 

k 


Wisi 

(A-l) 

Vi 

(A-2) 

If  one  defines  weight  vector  W  ■  [w,  .  .  .  wk]T  and  noise  vector 
N  -  [n^  .  .  ,  n..]T  (Ref  1:  590)  along  with  the  signal  vector  defined 

in  Section  I,  he  can  simplify  the  expressions: 


aWTS 

(A-3) 

wtn 

(A-4) 

The  next  step  in  determining  the  signal-to-noise  ratio  is  computing 
the  noise  power  in  the  output.  This  is  given  by  (Ref  1:  586  ): 

P„  ’  E[|v/] 

*  E  CvnV 

»  E  [(WTN)*(WTN)]  (A-5) 
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Noting  In  Eq  (A-5)  that  WTN  =  NTW,  one  can  further  write  (  Ref  1:586  ): 

Pn  =  (WT)*E[N*NT]W  (A-6) 

However,  E[N*N^]  =  M,  the  noise  covariance  matrix  (Ref  1:586  ).  This 
leads  to  the  following  equation  for  noise  power: 

Pn  =  (WT)*MW  (A-7) 

The  matrix  M  is  Hermitian  and  positive  definite,  so  it  can  be 
diagonalized  by  a  coordinate  transformation  (Ref  1:  586).  If  one  calls 
the  transformation  matrix  A,  the  signal  and  noise  vectors  in  the  trans¬ 


formed  coordinates  become: 

S  =  AS  (A-8) 

N  =  AN  (A-9) 

To  compute  Pn  under  this  transformation,  one  proceeds  as  follows: 

vn  =  WTN  (A- 10) 

Substituting  Eq  (A-9)  into  Eq  (A-10)  results  in: 

v  =  WTAN 

n 

=  (ATW)TN  (A-ll) 

Comparing  Eq  (A-ll)  with  Eq  (A-4)  implies  (Ref  1:  587  ): 

W  =  ATW  (A-12) 


Employing  the  transformed  vectors  gives  noise  power  as  (Ref  1:  587  ): 


(A— 13) 


Pn  =  (W')*E[N*N']W 

But  matrix  A  decorrelates  the  noise  components  and  equalizes  the  power  of 
each  (Ref  1:  587  ).  The  M  can  be  written  (Ref  1:  587  ): 

M  =  E[N*NT]  =  PnIk  (A-14) 

Choose  Pnsl  and  M  *  1^.  Thus,  Eq  (A-13)  becomes  (Ref  1:  587  ): 

Pn  =  (WT)*W  (A-15) 

Substituting  Eq  (A-12)  into  Eq  (A-7)  gives: 

Pn  =  (WT)*A*MATW  (A-16) 

Comparing  with  Eq  (A-15)  indicates  (Ref  1:  587  ): 

A*MAT  =  Ik  (A-17) 

M  =  (AV)-1  (A- 18) 

To  verify  the  control  law,  one  starts  with  the  optimum  weight  vector 
for  a  system  with  equal  and  uncorrelated  noise  components  (Ref  1:  587  ): 

WQ  =  yS*  (A-19) 

where  v  is  arbitrary. 

Under  the  transformation,  v$  is  given  by  (P.ef  1:  587  ): 
vs  =  a  WTS  (A-20) 

Applying  the  Schwarz  Inequality  to  Eq  (A-20)  yields: 

|vJ2i  l»|2[(viT)*W][(ST)*S] 


(S-21) 


Substituting  Eq  (A-15)  into  Eq  (A-21)  shows  that: 

/  <  |c|2[(ST)*S]  (A— 22) 

n 

This  determines  an  upper  bound  for  signal -to-noise  ratio.  If  one  shows 
that  Eq  (A-19)  produces  this  upper  bound,  then  he  demonstrates  the  Eq 
(A-19)  is  optimal  (Ref  1:  587  ). 

Putting  Eq  (A-19)  into  Eq  (A-20)  enables  one  to  compute  the  signal 
power  out  Jtef  1 :  587  ) : 


=  |ap(ST)*S|2 

-  |c|2|v|2[(ST)*S]2  (A-23) 

A  A 

The  noise  power  comes  from  replacing  W  in  Eq  (A-15)  with  WQ  in  Eq  (A-19) 
(Ref  1:  587  ): 

Pn  =  (WT)*W 

=  u*(ST)*yS 

=  |u|2[(ST)*S]  (A-24) 

Dividing  Eq  (A-23)  by  Eq  (A-24)  leads  to: 

p 

IT  *  M2[(ST)*S]  (A-25) 

n 

Consequently,  Eq  (A-19)  is  optimal.  However,  the  untransformed  array 
weights  are  needed,  not  the  transformed  weights  W. 


Appendix  B:  Second  Derivative  Approximation 


The  error  in  the  amplitude  distribution  revealed  by  the  spherical 
reflector  example  in  Section  III  arises  directly  from  inaccuracies  in  the 
approximate  computation  of  the  second  derivative  at  the  impact  point  on 
the  reflector  surface.  The  errors  in  amplitude  are  severe  because  the 
radius  of  curvature  of  the  reflector  depends  on  the  second  derivative. 

The  equation  for  the  surface  of  Figure  3-10  is: 

f (x)  =  Vr2  -  x2  -  h  (B-l) 

v 


where  R  *  200, 
•f“(x)  . 


h  *  150.  The  second  derivative  is: 
R2 

‘  (R2  -  x2)3/2 


(B-2) 


On  -50  <  x  <  50,  the  second  derivative  is  bounded  by  -0.00551  <  f"(x)  < 
-0.005;  so  f"(x)  is  always  small  in  magnitude.  Since  f"(x)  is  in  the 
denominator  of  the  expression  for  the  radius  of  curvature,  small  errors 
in  f"  can  generate  wild  variations  in  rc.  The  results  in  significant 
error  in  the  wavefront  radius  of  curvature  calculated  from  Eq  (3-7). 

This  is  the  root  source  of  the  amplitude  error. 

Unfortunately,  as  stated  on  page  205  of  Reference  3, "the  process  of 
differentiating  is  basically  an  unstable  process."  The  following 
formulas  were  used  (Ref  3:  206  ): 


f 


If 

0 


+  terms  of  order  h 


2 


(B-3a) 


applied  twice 


(B-3c) 


f:  * 


fi  -f-i 


Eqs  (B-3a)  and  (B-3b)  yielded  meaningless  results,  while  B-3c  produced 
data  that  clustered  around  the  distribution  found  from  using  Eq  (B-2)  to 
compute  the  second  derivative. 

Experimenting  with  double  precision  arithmetic  might  generate  better 

results;  however,  this  will  probably  not  solve  the  problem  in  general.  In 

the  most  general  case,  the  reflector  will  be  the  result  of  some  other 

numerical  algorithm.  Employing  double  precision  will  most  likely  yield 

only  a  long  string  of  meaningless  digits,  since  one  cannot  know  the 

reflector  surface  to  arbitrary  precision.  A  simple  experiment  conducted 

in  the  course  of  this  work  showed  no  difference  between  single  and  double 

precision  amplitude  results  for  the  10-wavelength  spherical  reflector.  In 

this  experiment,  the  second  derivative  was  calculated  with  double  precision, 

but  the  numerator  differences  in  the  difference  quotients  of  Eq  (B-3c) 

resulted  from  single  precision  arithmetic.  This  reflected  that,  in  general, 

one  cannot  know  f^  and  f^  to  arbitrary  precision.  Since  no  improvement 

followed  from  this,  one  must  conclude  that  errors  in  the  numerator  difference, 

2 

magnified  through  division  by  the  small  h  ,  produce  the  error.  Only  if 
algorithms  which  generate  the  reflectpr  grid  can  supply  very  great  precision, 
will  double  precision  arithmetic  prove  useful. 


Appendix  C:  Diffracted  Rays 


The  edge  diffracted  field  is  given  by  (Ref  12:  472): 

Ed(s)  =  -  D  Ei  A(s)  exp(-j^s)  (C-l) 

where 

s  ■  distance  from  diffraction  point  to  observation  point. 

E.-=  incident  field  phasor  at  diffraction  point. 

A(s)=  1/  s  (Spatial  Attenuation  Factor). 

D  =  diffraction  coefficient. 

The  s  is  computed  by  simple  analytic  geometry,  since  the  points  of  diff¬ 
raction  and  observation  are  both  known.  For  isotropic  feed,  the 
exp(-j^r1)/-^r^  where  r^  is  the  distance  from  feed  to  diffraction 
point. 

Calculation  of  D  is  somewhat  more  involved.  Figure  C-l  shows  the 
geometry.  Naturally,  the  diffracted  ray  in  question  is  the  one  which 
intersects  the  aperture  point  under  consideration.  The  angle  0  must  be 
computed  based  on  the  location  of  the  aperture  point  with  respect  to  the 
point  of  diffraction.  The  0',  however,  depends  only  on  the  feed  location 
and  the  edge  location.  The  angles  are  measured  with  respect  to  the 
tangent  line  to  the  curve  at  the  edge.  Once  the  angles  have  been  determ 
ined,  one  should  compute  (Ref  12:  473-475): 

S  rl 

L  3  - L  (C-2a) 

s  +  rx 

N+  .  7T+  (0  +  0') 

N+  "  -  rounded  (C-2b) 

47 T 

N+  =  7T  +  (0  -  0') 


rounded 


(C-2c) 


rounded 


(C-2d) 


TT  +  (0  +  0* ) 


-  7T  +  (0  -  0') 

N_  - - rounded 

4 


( C-2e ) 


The  N's  are  integers  found  by  rounding  the  given  expressions  to  the 
nearest  integer  (Ref  12:  475).  They  are  in  turn  used  to  compute  (Ref  12 


a+(0  +0')  =  2  cos2 


477 Nj  -  (0  +  0') 


(C-3a) 


a"(0  +  0' )  =  2  cos' 


4TT  N+  -  (0  +  0') 


4  .  2  r  4?TN+  -  (0  -  0')] 

a+(0  -  0')  *  2  cos^  - I - 


(C-3b) 


(C-3c) 


Defining  the  following  is  helpful  (Ref  12:  474) 


Q  =  -exp(-jTT/4)/(4-^2Tr  S  ) 


[t] 


Also,  F(X)  as  follows  is  needed  (Ref  12:  474): 

,00 

F(  X)  =  J2  |V^”|  exp(  jX)  lex  p(-jt2)  dt 


(C-5) 


With  all  the  preceeding  definitions,  one  can  calculate  the  diffraction 
coefficient  (Ref  12:  474): 

D  =  Q  [  Kx  F(f  L  a+(0  -  0*))  +  K2  F(  jk  a"(0  -  01)) 

-  K3  F(f  L  a+(0  +  0'))  +  K4  F(^L  a“(0  +  0'))]  (C-6) 


The  reader  should  note  that  Eq  (C-6)  applies  only  to  the  problem  at 
hand,  a  straight  edge  with  an  incident  field  parallel  to  the  edge.  The 
reader  should  consult  Reference  12  and  available  literature  for  any 
other  problems. 


Diffracted  Ray 


Tangent  Line  to  Reflector  Surface  at  Edge 


**  / 


Appendix  D»  Programs 


This  appendix  contains  programs  used  in  this  work.  The  language  is 
77  standard  Fortran.  The  purpose  of  each  program  or  subroutine  is  noted 
in  the  listings  themselves. 
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* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


Sub  routine  AFIELD  <  RZ  *  NGRID .  r I  ,  YA ,  i-  RAYS  ,  Et ,  YEA ,  Nr  > 

APTFLD  is  the  controlling  routine  For  the  trac:r  : 

method  and  will  calculate  an  approximate  aperture  field. 

RZ  is  the  grid  of  reflector  points,  given  as  distance  above 
the  aperture  plane* 

NGRID  is  the  number  of  grid  points* 

TI  and  YA  are  working  arrays.  The  NRAYS  should  be  chosen  as 
some  appreciable  fraction  of  NGRID,  say  1/4. 

EA  is  the  array  of  approximate  aperture  field  values  an a 
YEA  is  the  array  of  aperture  points.  The  NAP  is  the  number 
of  aperture  points. 

Inputs*  RZ, NGRID, NRAY«,NAF 


Real  RZ<NGHIB> ,TI(NRAYS> ,YA(NRAYS,6) , YEA (NAP) 

Complex  EA(NAP) 

Common  /BOUNDS/  RF YM IN, RFYMAX , YEAMIN , YE AM AX 

Common  /CNSTNT/  PI, BET A 

HG  »  <RFYMAX  -  RFYMIN) /FLOAT (NGRID-l > 

HA  =  <YEAMAX  -  YEAMIN) /FLOAT (NAP-1 ) 

Do  1  I»1,NAP 

YEA< I )  »  YEAMIN  +  HA#FLQAT ( 1-1 ) 

EA<  I )  »  <0.0,  0.0) 

1  Continue 

THIN  »  ATAN< RZ< NGRID )/RFYMAX) 

TMAX  »  PI  -  ATAN<ABS<RZ<  D/RFYMIN)) 

DT  »  <TMAX  -  TMIN)/FL0AT<NRAYS-1 ) 

Dp  2  I»l, NRAYS 

TI < I )  a  TMIN  +  DT*FL0AT<I-1> 

Call  TRACE <HG, RFYMIN, RZ,NGR ID, TI < I ) , 

1  YA<I,l),YA<I,2),YA<I,3),YA<I,4>,YA<I,w),YA<I,A); 

2  Continue 

Do  7  1=1, NAP 

Q  «  YEA< I >  -  YA< 1,1) 

If <ABS<Q) »LT» 1 .OE-4) then 

Call  FIELD <YA< 1,1),YA<1,2) ,YA< 1 ,3) , YA< 1 ,4) , YA< 1 ,5) , YA< 1 ,6) , 

1  EA< I ) ) 

Endif 

Do  6  K»2, NRAYS 

QQ  a  YEA<I )  -  YA(K, 1 ) 

If < ABS<QQ) .LT.l. OE-4) then 

Call  FIELD < YA<K, 1 ) , YA<K,2) ,YA<K,3) ,YA<K,4) ,YA(K,5> , YACK, A) 
1  EA< I ) ) 

Go  to  5 
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w  «r  in  -o  rs 


If  <  <ABS(Q>  .LT.l  .OE-4)  .OR.  (ABS<CC5)  .  J 
If  (QHcQQ  ♦  6T  ♦  0 . 0  )  Go  to  5 
TLOW  *  TI  <  !<-l ) 

THIH  =  TT (K> 

01  =  a 

GG1  =  QG 


—4  '  !>J  to  3 


Dec  1  00*01  1983  aperture  Page  2 


Do  3  JJ=1,100 

TM  »  <TL0W  +  THIH)/2*0 

Call  TRACE ( HG ? RFYMIN ? R2 ? NGR1B , TM , YAPT , YT , Z I , THR , DZDY , D2ZDY2 ! 
QQQ  *  YEA( I >  -  YAPT 
If (ABS(QGQ)  *LT *  1.0E-4)  Go  to  4 
If<Gl*GGQ  »GT «  0.0)then 
TLOW  =  TM 
G1  =  GGQ 
Endif 

If (GGG&GG1  .GT*  0.0>then 
THIH  -  TM 
QQ1  ■  QQQ 
Endif 
Continue 

Call  FIELD (YAPT ?YI ?ZI * THR? DZDY ?D2ZDY2  ? EA< I ) ) 

Q  *  00 
Continue 
Continue 
Return 
End 

Subroutine  TRACE ( H ? YMIN * Z ? NG * THETA , YAPT ? Y I ? ZI ? THR , DZDY , D2ZDY2 ) 

TRACE  does  the  actual  ray  tracing  for  given  incident  ray 
angle  with  the  horizontal* 

H  is  the  reflector  grid  spacing  and  YMIN  is  the  lowest 
point  on  the  aperture* 

Z  is  the  reflector  grid  with  size  NG* 

THETA  is  the  angle  that  the  incident  ray  eakes  with  the 
horizontal?  Measured  counter-clockwise  positive. 

YAPT  is  the  intersection  point  of  the  reflected  ray  with 
the  aperture  plane* 

YI  and  ZI  specify  the  impact  point  of  the  incident  ray. 

THR?  DZDY?  D2ZDY2  are  the  angle  of  reflection  with  respect 


•e*  .  *  *  V V  *  *  ■  ' 


*  *  *  * 


Hr  to  the  surface  normal,  the  first  ienvol-.ve,  an n  the  second 

derivative  respectively ♦ 


Inputs:  HrYMTNrZ»NOf THETA 

Real  Z<NG).K- 
Common  /C^GTr-T/  PI. BETA 
If <ABS< THETA  -  PI/2.0)  »LT .  l.CE-4)then 
YT  »  0.0 

ZI  *  Z< <NG-l)/2+l> 

Call  DERV<<NG-l)/2  +  l,Z,NG,W,£iZBY,D2ZDY2) 
Call  BOUNCE <  Y I , ZI , DZDY , THETA , THR ,  Y  APT ) 
Return 
End  if 

K  *  TAN (THETA) 

YNH1  -  YHIN 

GNM1  -  K*YNIN  -  Z<1> 
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If ( ABS<GNM1 )  .LT.  1.0E-4)  then 
YI  «  YMIN 
ZI  »  Z(l) 

Call  DERV(1,Z,NG,H,DZDY,D2ZDY2> 

Call  BOUNCE <  YI , ZI , DZDY , THETA , THR , YAPT ) 
Return 
Endif 

Do  3  N«2,NG 
YN  -  YNM1  +  H 
GN  »  K*YN  -  Z(N) 

If < ABS(GN)  .LT.  1.0E-4)  then 
YI  *  YN 
ZI  -  Z(N) 

Call  DERU(N»ZfNG»Hf DZDYrD2ZDY2> 

Call  BOUNCE < YI , ZI , DZDY f THETA , THR , YAPT ) 
Return 
Endif 

If <GNM1*GN  »GT»  0.0)  Go  to  4 
F  *  GN/(GN  -  GNM1 ) 

YI  *  YN*( 1 .0  -  F)  +  F*YNM1 
ZI  •  Z(N)*< 1 .0  -  F)  f  F*Z(N-1) 

Call  DERV < N- 1 » Z » NG , H , DZNH 1 , D2ZNM 1 ) 

Call  DERV(NvZfNGfHf DZDYNf D2ZN) 

DY  -  YI  -  YNM1 

DZDY  -  DZNM1  4  D2ZNM1*DY 

D2ZDY2  «  D2ZNM1  +  (D2ZN  -  D2ZNM1)*DY/H 

Call  BOUNCE < Y I . ZI , DZDY , THETA , THR , Y APT  > 

Return 

4  YNH1  ■  YN 


9* 


WO.iUllL 


rrrrrr*  T  7T 


GNM1  a  GN 
t>  Continue 
E  nd 


TTSTrTfT^TTry^ 


* 

»JY 


Subroutine  BOUNCE  <  Y I , ZI .DZDY .THETA . 

BO UMr,£  serve:?  TRACI- .  T!-e  Yl.Z 


HR  f  Y APT  > 

T  *  T  I-  •  ;>  •  J  .”?'f 


Common  /CN3TNT/  PI ? BETA 
If  (ABSdZDY)  . LT ♦  i.OE-6)  then 
PHI  a  PI/2.0 
Go  to  1 
End  if 

If (DZDY  .LT.  0.0)  then 
PHI  ■  ATAN( 1 »0/ABS(DZDY) ) 

Go  to  1 
Endif 

PHI  ■  PI  -  AT  AN ( 1 » 0/DZDY  > 

1  PSI  -  2.0*PHI  -  THETA 
THR  «  ABS(PHI  -  THETA) 

If < ABS(PSI  -  PI/2.0)  .LT.  1.0E-4)  then 
YAPT  »  YI 
Return 
Endif 

If (PSI  .LT.  PI/2.0)  then 
YAPT  »  YI  -  ZIZTAN(PSI) 

Return 


.  .V: 
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Endif 
YAPT  -  YI 
Return 
End 


+  ZI/TAN(PI-PSI) 


Sub  routine  FIELD  <  YAPT ? Y I ? ZI ? THR , DZDY ? D2ZDY2 ? EAPER > 

FIELD  calculates  the  contribution  of  the  reflected  ray 
to  the  aperture  field. 

YAPT? YI ?ZI ?THR?DZDY?D2ZDY2  are  inputs.  They  are  the 
intersection  of  reflected  ray  with  aperture?  impact 
point?  angle  of  reflection  with  respect  to  surface 
normal  at  impact  point?  and  first  and  second  derivatives 
at  impact  point  respectively. 

Complex  EAPER 
Real  L 

Common  /CN8TNT/  PI ? BETA 
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BY  *  YAP7  -  Y I 

RHO  -  SORT  ( YIIcYI  +  21*21) 

RHOA  *  SORT <21*21  +  DY*DY> 

L  ■  RHO  +  RHOA 
PHI  »  -RETAIL 
A  *  1 . 0/SQR7 < SHG) 

I f  <  ABS  '  B22DY2  )  .  LT .  1.  OE-6 )  then 
A  =  A*SnR7<RHC/CRH0  +  RHOA ) ) 

Go  to  1 
Endif 

RC  -  <1*0  +  DZDYKDZDY ) ##1 . 5/D2ZDY2 
U  «  AB3< 1 . O/RHO  +  2.0/(RC*C0S(THR>)> 

If(U  .LT •  1 *0E-5)Go  to  1 
RCW  -  1.0/W 

If (D2ZDY2  .GT.  0.0)  A  ■  A*SQRT(RCW/<RCW  +  RHOA)) 

If (D2ZDY2  *LT.  0.0)  then 

If ( ABS(RCU  -  RHOA)  .LT.  1.0E-5)Go  to  1 
If (RHOA  .LT.  RCW)  A  *  A*SQRT(RCW/(RCW  -  RHOA)) 

If (RHOA  .GT.  RCW)  then 

A  »  A*SQRT( RCW/ (RHOA  -  RCW>> 

PHI  *  PHI  +  PI/2.0 
Endif 
Endif 

1  EAPER  =  EAPER  ¥  A#CEXP(CMP!.X(0.0,PHI  >  > 

Return 

End 

Sub  routine  DERV  <  N , Z , NG , H , DZDY , B2ZDY2 > 

* 

*  DERV  simply  approximates  the  derivatives  of  the  reflector 

*  grid  at  grid  points  using  established  numerical 

*  differentiation  formulas. 

* 

*  N  is  the  point  in  questionf  while  Z  is  the  reflector  grid 

*  NG  is  its  dimension.  H  is  the  spacing.  These  are  all  input. 

* 

Real  Z(NG) 
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If (N  *EQ.  1)  then 

DZDY  -  FR2(Z(N),Z(N+l),Z(N+2)fH> 
DZDY1  -  C2(Z(N) » Z(N+2> »H) 

DZDY2  «  C2(Z(N+1) »Z(N+3) »H) 

D2ZDY2  ■  FR2(DZDY» DZDY1 »DZDY2f H) 
Return 
Endif 

If (N  .EQ.  NG)  then 

DZDY  »  FR2(-Z(N),-Z(N-1>,-Z(N-2>,H> 


* 


mmwrrzwrrz  *1* £.»■!*!  ■  r  "  rf?*”.*  ■:  i  v  ^  ■  i1  v  n**  ^ ^ vy T» ^ •.-■  •.■« t*  ■>  .■»  o V" 

UZBYM1  *  C2  <  Z  <  N-  2 )  *  2  >  N )  *  H ) 

B2BYM2  -  C2<Z<N-3),Z(N-l)fH) 

B2ZBY2  »  FR2 < -B2BY , -BZBYM 1 , -DZDYM2 , M ) 

Return 
End  if 

If  ( <N  *EQ»  2 >  . CR .  ( N  . E2 .  N3-  1  > th »• 

02BY  a  22  ;  Z  ( N- _  ;  -  Z  <  w-  :• )  «p») 

If  ( N  .t:Q.  2)  than 

DZDYHl  *  FR2iZ(N-l).Z(N>fZfN+l)»H) 

BZBY1  *  C2  <  Z ( N  > 1 2  <  N-r2  >  » H ) 

D2ZDY2  »  C2(BZBYM1,BZBY1,H> 

Else  If (N  *EQ.  NG-1 >then 
DZDYM1  *  C2<Z(N-2> ,Z(N> ,H> 

DZDY1  *  FR2<-Z<N+1),-Z<N>,-Z<N-1),H> 

D2ZDY2  ■  C2<DZDYMl,DZDYlrH) 

End  if 
Return 
End  if 

DZDY  »  C2<Z(N-l)fZ<N+l)»H) 

BZDYM1  =  C2(Z(N-2> ,Z(N> ,H> 

DZDY1  -  C2(Z(N> ,Z(N+2> ,H>  j 

D2ZBY2  *  C2(BZBYM1,DZBY1,H> 

Return  I 

End 

Function  FR2<X0,X1 ,X2,H) 

FR2  *  <-3.0#X0  +  4.0*X1  -  X2>/(2.0*H> 

Return 

End 

Function  C2(XM1,X1,H) 

C2  a  (XI  -  XM1>/<2.0*H> 

Return 

End 
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* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 
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Sub  routine  Ap7INT  i  ~  ,  F .  X  ,  "  --lIT A .  N» 


IN  *  YMAX . 


APTINT  computer  the  far  field  for  a  g  -  ven  aperture 

distribution  ♦ 

E  and  P  are  the  amplitude  and  phase  respectively  of  the 
aperture  distribution*  Input 

X  is  a  complex  scratch  array* 

NA  is  the  dimension  of  the  E»P,X. 

F  is  the  pattern  that  results* 

THETA  is  the  angle  from  boresight  associated  with  the 
corresponding  F  value* 

YMIN  and  YMAX  are  the  dimensions  of  the  aperture* 

TMIN  and  TMAX  give  the  range  of  angle  for  which  the  potter, 
will  be  calculated* 

Real  E<NA>  ,P<NA),F(NP>, THETA < NP ) ,MAX 

Complex  X<NA) f Q 

Common  /PIBETA/  PI, BETA 

H  »  < YMAX  -  YMIN)/FL0AT<NA-1> 

DTHETA  *  <TMAX  -  TMIN>/FL0AT<NP-1> 

Do  2  I-1,NP 

THETA < I )  -  TMIN  +  DTHETA*FLOAT ( 1-1 > 

STH  -  SIN(THETA< I )  ) 

Do  1  J*l,NA 

Y  »  YMIN  +  H*FL0AT<J-1> 

X(J)  -  E<J>*CEXP<CMPLX<0.0,BETA*Y*STH  +  P<J>)> 

1  Continue 

Call  SRL<X,NA, YMIN, YMAX, Q> 

F< I )  *  CABS ( Q ) 

2  Continue 
MAX  *  F< 1 ) 

Do  3  1*1, NP 

If <F< I >  *QT*  MAX) MAX  *  F<I) 

3  Continue 

Do  4  1*1 ,NP 

F( I )  *  F < I ) /MAX 

THETA Cl)  ■  THETA< I )<< 180* 0/P I > 

4  Continue 
Return 
End 


Subroutine  SRL<F,N,XMIN,XMAX,Q) 

* 


J  M"J  ^  9  FJlT.  &.'  r:v:  ^  •"-  *Vr-„» ■•■ : .■  r 


•SSMSSRl 

*  Simpson's  Rule  integration  routine  asr-oc  iated  with  Al*i  IN’*' » 

* 

Complex  F(N),Q 
H  *  (XMAX  -  XMIN) /FLOAT (N-i ) 

!.  -  M  -  2  ' 

TF<  <  <N~1)/2)S2  ♦  NF  »  <  N— 1 )  )'_  =  L  -  3 
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Q  »  <0,0, 0.0) 

Do  1  1=1, L, 2 

Q  -  Q  +  <F< I )  +  4 .  0#F ( 1  +  1 )  +  F < 1+2) ) 

1  Continue 
Q  =  H*G/3,0 
If (L  ,EQ,  N-2)Return 

Q  *  Q  +  3,0*H#<l-<N-3)+3,0#F(N"2)+3,0#F < N-l )  +F <N>  >  /8 .0 

Return 

End 


0 
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Program  CNURT 

* 

*  This  program  converts  the  rho  versus  osi.  da-. a  f  ~c-r.  .n 

beamshaping  formulas  into  a  equi-spacea  x-y  format. 

Real  RH0< — dimension  of  rho — )»PSI< — dimension  of  psi — ) 
0pen<l,File*'name  of  input  file') 

Open (2, File* 'name  of  output  file') 

NP  *  — number  of  input  dato  points — 

NG  ■  --number  of  output  data  points — 

YMIN  »  — lower  limit  of  aperture — 

YMAX  *  — upper  limit  of  aperture — 

HG  *  (YMAX  -  YMIN) /FLOAT (NG-l > 

Rewind  1 
Do  100  K=1,NP 

Read ( 1 f 5)  PSI <K) ,RH0<K) 

100  Continue 

Do  6  1=1,  NR 

Y  =  YMIN  +  HG#FL0AT(I-1) 

Q  *  Y  -  RH0< 1 ) tSIN <  PSI < 1 ) ) 

If ( ABS  <  Q )  .LT »  1.0E-5)then 
P  *  PSI<1) 

R  *  RH0< 1 ) 

Go  to  4 
Endif 

Do  3  J=2»NP 

QQ  *  Y  -  RHO(J)*SINtPSICJ>) 

If < ABS(QQ)  .LT.  1 .OE-5) then 
P  *  PSI<J) 

R  *  RHO(J) 

Go  to  4 
Endif 

If(Q*QQ  .GT»  0.0) Go  to  2 

DRDP  *  <RH0< J)  -  RHO(J-l) )/(PSI( J)  -  PSI(J-l)) 

01  *  Q 
Q2  *  QQ 
PI  »  PSI(J-l) 

P2  *  PSI(J) 

1  P  *  <P1  +  P2)/2.0 

R  *  RH0< J-l )  +  DRDP*(P  -  PSI(J-l)) 

X  »  Y  -  R*SIN(P) 

If <ABS(X)  .LT.  1.0E-5)Go  to  4 
If (Q1*X  .GT.  0.0)then 
Q1  -  X 
PI  *  P 
Go  to  1 
Endif 
Q2  *  X 
P2  -  P 
Go  to  1 


fc.lA.ia.i  a  <k  fcV  t  .  y*  «-  V  .»  ’j  ' 


:  i  ro 


Q  •*  QQ 
Con  tin ue 
Z  3  3*CQS<P> 
Write <2 >5)  Y,Z 
5  Fori.iat<2Fl2»5) 
f:  Continue 
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Stop 

End 
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